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We analyze the potential dark matter implications of LHC events with missing transverse mo¬ 
mentum and a resonance, such as a Z', decaying to a pair of jets or leptons. This final state contains 
significant discovery potential, yet has not yet been examined in detail by the LHC experiments. We 
introduce models of Z' production in association with dark matter particles, propose reconstruction 
and selection strategies, and estimate the sensitivity of the current LHC dataset. 

PACS numbers: 


INTRODUCTION 

As the Large Hadron Collider (LHC) resumes opera¬ 
tions this year after a major upgrade and a half-decade 
of data taking, a central area of focus will be the search 
for physics beyond the standard model (SM). While the 
LHC will have sensitivity to many models inspired by 
theoretical extensions or generalizations of the SM, the 
search for dark matter is of particular interest due to 
the well-established fact of its existence JJ. The collider 
detection of dark matter is a cornerstone of the effort 
to elucidate and obtain evidence for the particle nature 
of dark matter, and is complementary to astrophysical 
methods of detection. 

Searches for dark matter production at the LHC rely 
on the production of a visible object X recoiling against 
the missing transverse momentum ($t) from the invisible 
dark matter particles. Cases where X is a SM particle 
such as g/q pH?], W [SHE], Z [El E5] , H [HH5] , 7 (3] 
HE], or a heavy quark HMD have been considered. 
For a review of simplified models for dark matter at the 
LHC, see Refs. [221125] . 

In this paper, we present a new mechanism for dark 
matter production at the LHC, where the visible object 
is itself a new particle, a Z' boson. We propose examples 
of models giving rise to a signal of Z' + Ifi t, where the 
Z' boson can decay to pairs of charged leptons {I + I~) or 
to pairs of quarks leading to jets (jj), and is therefore 
distinguishable as a resonance in the dilepton or dijet 
mass spectrum. In each case, we study the sensitivity 
of the LHC in this channel, and compare with existing 
searches for the Z' without a requirement of large 

The models here specifically target the production of a 
new Z' which is present in a hypothetical, non-minimal 
dark sector. New Z' bosons arise in many extensions to 
the SM [23], and the possibility of dark matter coupled 
to a Z' has been explored extensively in the literature, 
including in the context of the LHC (see, e.g. [3] [251 - 


[55]). It should be noted that the experimental signature 
of a dijet or dilepton resonance plus missing transverse 
momentum does not require a Z' : other possibilities, in¬ 
cluding new scalar resonances or colored resonances, are 
natural directions to explore. 

In addition to extending the current program of X +^t 
studies, the models presented here point to final states 
whose LHC data remains unexamined and which are nat¬ 
ural generalizations [551 of previously performed searches 
for Z/W + with Z —► U or Z/W — > jj. These data 
therefore contain real, untapped discovery potential, in¬ 
dependent of theoretical interest in models of dark matter 
involving Z' bosons. 

The models considered here are also examples of dark 
sector signals that, to some extent, could be hidden in 
existing ^ T -based searches. Searching specifically for a 
dijet or dilepton resonance reduces the backgrounds and 
could give a strong hint of new physics. Furthermore, 
many searches have been optimized for new high-mass 
particles. For the examples below, we find that the most 
unconstrained parameter space is for relatively light Z' 
states, those with vriz 1 below 100 GeV, where current 
LHC searches have low efficiency. 

In the following, we first review experimental con¬ 
straints on Z' gauge bosons and then describe several 
models of Z' + dark matter production. The range of 
Z' mass explored is 50-800 GeV: for lower masses, dijet 
masses would be more difficult to reconstruct due to a 
smaller angular separation in the partons. Work on LHC 
signals of Z' + with lower values of mz' will appear 
elsewhere [35] , while related work focusing on leptonic Z' 
decays plus missing transverse momentum can be found 
in Ref. [ 551 . 

We consider two models with a minimal set of renor- 
malizable interactions: dark-Higgsstrahlung from a Z ', 
with the dark Higgs decaying invisibly; and a dark sector 
with two states Xi ,2 that couple off-diagonally to the Z'. 
We also study the case where the production of the dark- 
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sector states is through a higher-dimension operator. We 
analyze the sensitivity of the current LHC run to these 
models in jj + $t and It + $t final states, and compare 
to existing constraints. For the renormalizable models, 
the Z' + search has better sensitivity than direct res¬ 
onance searches only for low Z' masses. In the operator 
case, it is possible to probe the scale of new interactions 
to around a few TeV. 

CURRENT CONSTRAINTS ON Z' BOSONS 

For simplicity, we assume a U(l)' where the Z' has 
universal vector couplings to SM quarks: 

CD-' 52 g q <rfqZ'>*. ( 1 ) 

i 

The couplings above are the same as for gauged baryon 
number U(1)b with g q = gs/S, where anomaly cancel¬ 
lation could be achieved with additional heavy quarks 
or with chiral matter in a dark sector. This possibility 
has been studied in detail in the context of dark mat¬ 
ter (e.g., [37, [HE]). However, we do not assume that the 
gauged baryon number is the origin of the U(l)'. For ex¬ 
ample, it is possible that the Z' couplings to SM fermions 
are generated by higher-dimensional operators [3Hj while 
the dark sector states are directly charged under the Z'. 

When we consider dilepton searches, we will introduce 
additional free parameters for couplings of the Z' to lep¬ 
tons. Since the production of the Z' does not depend on 
the lepton coupling (except through the dependence on 
the Z' width, which we neglect) we present constraints 
from dilepton resonances searches simply in terms of the 
Z' branching ratio to the appropriate lepton flavor. 

Although we will not impose any relationship between 
the Z' coupling to quarks, leptons, or dark sector parti¬ 
cles, one natural possibility is that of kinetic mixing m, 
where a mixing of Z' and hypercharge generates cou¬ 
plings of the Z' to SM fermions. Since the natural size 
of the couplings is small in this case ( 10~ 2 or less), we do 
not consider this for the models that rely on dark matter 
production via the Z' couplings to quarks. However, this 
gives a simple way for the Z' to decay to visible states 
in our last model, where the Z' is only produced in the 
decay of dark sector states. 

The range of Z' mass explored here is 50-800 GeV. For 
heavier masses, constraints from dilepton or dijet reso¬ 
nance searches are precisely where LHC searches excel 
since backgrounds are relatively low. The Z' + sig¬ 
nature has additional particles produced along with the 
Z' and so has a smaller rate than direct Z' production; 
therefore we expect it to be a less sensitive probe of the 
models in the high mass regime. Meanwhile, a low mass 
Z' decaying to quarks is difficult to resolve as separate 
jets; however, it is possible that this regime could be stud¬ 
ied by employing jet substructure techniques [351 W- 


Dijet Constraints 

Direct dijet resonance searches constrain a Z' coupling 
to quarks. We take limits on g q as a function of Mz< from 
Ref. [32]) which compiles experimental results down to 
Mz' = 140 GeV. Here, the lowest mass region was cov¬ 
ered by UA2 [43] with integrated luminosity of 10.9 pb -1 . 
At lower Z' mass, dijet resonances are more difficult to 
constrain due to the large QCD background. Data on 
the dijet spectrum down to rrijj = 48 GeV have been 
published by UA2 [33] (4.7 pb -1 ) and down to 60 GeV 
from CDF [45] (26 nb _1 ). While a reanalysis of the data 
would be needed to obtain limits on new resonances, we 
estimate that the UA2 dijet limits continue to weaken 
below 140 GeV, reaching g q < 1 at Mz> = 50 GeV (see 
also [IS]). 

Future LHC analyses may be able to provide more ro¬ 
bust coverage of the low mass Mz> region. This was stud¬ 
ied in Ref. |28], which considered associated Z' searches, 
such as a Z' in addition to a Z, 7 , or jet. Using the 
additional object in the final state to improve trigger ef¬ 
ficiency, it was shown that LHC searches can have sensi¬ 
tivity even down to Mz> ~ 50 GeV and couplings com¬ 
parable to or better than the estimated UA2 dijet limits. 


Dilepton Constraints 

A Z' coupling to electrons is strongly constrained by 
LEP measurements m- In the first two of our models, 
we will focus on the possibility that the Z' has suppressed 
couplings to electrons but 0 ( 1 ) branching ratio of the 
Z' to muons (for example, see Ref. [3H] and references 
therein). 

Then if the Z' has a preferred coupling to muons and 
to quarks, a hadron collider can give interesting limits 
relative to the LEP precision measurements. As a direct 
comparison to Z' + searches, we consider constraints 
from searches for dimuon resonances. Limits are available 
from the CDF collaboration [39] down to Mz< = 100 
GeV, while ATLAS [50] and CMS |5T] limits extend down 
to Mz' = 150 GeV and Mz> = 300 GeV, respectively 1 . 
Here published results are not available below Mz due 
to the large Drell-Yan background. 

Below the Z-pole, Ref. [53] showed that LHC mea¬ 
surements of the Drell-Yan spectrum at low invariant 
mass [55] can be used to set strong constraints on a Z' 
coupling to quarks and muons. The recast of the data 
leads to constraints on couplings at the 10 -3 —10“ 2 level. 


1 LHC searches for the SM Higgs decay to dimuons can also be re¬ 

cast to place constraints down to Mz> = 110 GeV [521153] : we do 
not consider these analyses as they are not directly applicable to 
our models, and would not qualitatively change our conclusions. 
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(In the context of kinetic mixing, the current constraint 
is e < 0.012 and can reach e = 5 x 10 -3 for a binned 8 
TeV LHC analysis at Mz> = 50 GeV.) 

We also consider small, universal couplings of the Z' 
to all of the charged leptons, as in the case of kinetic 
mixing. As discussed above, this will be most relevant 
in our third model (Inelastic EFT) where the Z' may be 
very weakly coupled to SM fermions. 


Other Limits 

A light Z' coupling to quarks contributes to the Z 
hadronic width through Z — » qqZ' —» 4 j and through 
a Zqq vertex correction [551EZ!. Applying the results of 
Ref. [SSj to the most recent measurement of Rz = F(Z 
hadrons) /T(Z —y g, + g~) = 20.785±0.033 [58], this places 
a constraint of g q < 0.4 — 0.6 for Mz> = 50 — 140 GeV 
where there are no dijet resonance constraints. 

Finally, even if a Z' couples only to quarks, kinetic 
mixing of the Z' with Z, 7 can be generated at one- 
loop. There are strong constraints on this kinetic mixing 
from precision electroweak measurements [57][HD. giving 
a bound e < 0.02 for Mz< -C Mz- Since the kinetic mix¬ 
ing parameters are model-dependent, we do not examine 
this constraint any further, except to note that it is par¬ 
ticularly strong for Mz> ~ Mz and so any model in this 
case would have to be particularly tuned. 


MODELS OF Z'+$t PRODUCTION 
Dark Higgs 

A model with a new Z' naturally comes with its own 
scalar (or set of scalars) responsible for spontaneous sym¬ 
metry breaking. Suppose there is a new massive scalar 
that couples to the Z', which we call the dark Higgs, 
ho. Analogous to the SM process of Higgs-boson radi¬ 
ation from a W or Z, the new scalar is radiated by the 
Z' in a dark-Higgsstrahlung process. If this new dark 
Higgs boson additionally couples with invisible states 2 * , 
its primary signature could be as shown in Fig. [T] 

As a minimal model for this process, we introduce a 
new U(l)' with a charged scalar field $£> and an invisible 


2 Another possibility is that the Z' decays to dark matter, while 

the dark Higgs decays to SM states through mixing with the SM 
Higgs. Then the monojet search channel would also be sensitive 

to the model. 



FIG. 1: Diagram of the production of a Z' in association 
with a dark Higgs boson (ho) which decays into two stable 
dark states, It is assumed that hn is lighter than 2A I Z ' 
and decays with 100% branching to the invisible states. 

singlet scalar <\>x'- 

a d\d^ d \ 2 + vI\$d\ 2 - \d\$d\ 4 - \(F 

+ 2 ^ d ^x) 2 - ^x\$d\ 2 4> 2 x - V{(j>x)- (2) 

The dark Higgs field 4>£> = -^{vd + ho) obtains a vev 
vd, giving mass to the Z'. The masses of the dark scalars 
hjj and 4>x are fixed by the scalar potentials, and the 
Z' couplings to quarks are as in Eq. [T] Furthermore, if 
m x > 100 GeV or is very close to m-h/ 2 , it is straight¬ 
forward for <f>x to be a good thermal relic dark matter 
candidate if a scalar Higgs portal coupling is added to the 
Lagrangian in Eq. [ 2 ] [60]. However, we do not require <px 
to be a thermal relic as this would impose a restriction 
on M ho . 

The coupling of hp with the new gauge boson is 

Q h g z M Z 'h D ZlZ'» = g hD Mz'h D Z' ll Z'^, ( 3 ) 

where Qh is the charge of $£), which is a free parameter 
that we absorb by defining the effective coupling gh D - 
The dark Higgs can decay dominantly to the invisible 
4>x states through the Xx coupling, which we can take 
to be 0(1). Meanwhile, decays of hp -A Z'Z'* will be 
suppressed as long as m.h D < 2 Mz>. We assume the 
mixing of hn with the SM Higgs is small. 

As discussed in the previous section, we will take the 
SM charges under the Z' to be a separate free parameter, 
in order to be as general as possible. In considering sig¬ 
natures with dijets plus missing transverse momentum, 
we consider only the coupling to quarks; for dilepton plus 
missing transverse momentum signals, we will focus on 
the possibility of a non-zero branching fraction to muons. 

The masses M hr> and M z > are independent quantities 
in the model, though they are set by the same scale vd, 
with Mh D /Mz ' = V^Xd/ gh D - Note that since the Z' + 
signal due to the process shown in Fig.[l]favors larger 
gh D , this implies that for perturbative couplings the dark 
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FIG. 2: Diagram of the production of X 1 X 2 , followed by 
decay of the heavier dark sector state X 2 to Z 1 + x 1 , where 
Xi is a possible dark matter candidate. 


Higgs cannot be much heaver than the Z '. In order to 
capture most of the effects of different particle masses, 
we simply consider here two benchmark scenarios. In 
the “light” Mh D benchmark case, we set: 


Mho 


M z > , M z > < 125 GeV 
125 GeV , M z , > 125 GeV, 


(4) 


In the “heavy” Mh D benchmark case, we set 3 : 


M ho 


125 GeV , M z , < 125 GeV 
M z , , M hD > 125 GeV. 


(5) 


Light Vector 

When the Z' is relatively light, it can be produced in 
the decays of dark sector states 4 . An example is given in 
Fig- ® where the Z' possesses off-diagonal couplings to 
dark sector states \i and Xi- If the mass splitting be¬ 
tween the two states is larger than M z < , the heavier state 
(X 2 ) can decay to an on-shell Z' and a Xi- Meanwhile xi 
is stable and a dark matter candidate. 

As a concrete example, we consider a Z' coupled to a 
new fermion which has both Dirac and Majorana masses. 
The fermion \ initially has a Dirac mass Md and vector 
coupling with respect to the Z'. A Majorana mass can 
be generated from the vev of a U{ 1)' Higgs through an 
interaction y x $> x xx C i so that 


C D X(}JJ> ~ M d )x - ^~(XX C + h.c.). (6) 


3 For the lowest mass point considered Mz> = 50 GeV, the decay of 
hf) —>■ Z'Z' is kinematically allowed; for simplicity we continue 
to fix the h ~£) invisible branching fraction to 1. 

4 Alternatively, the Z' can be produced as radiation from off-shell 

dark sector states [35] |36j. 


This will lead to two Majorana states x' 1,2 with masses 
Mi 2 = | M m ± Md\. The interaction with the Z' is off- 
diagonal and can be written as: 

Y Z h (x27 /i 7 5 Xi + Xi7 / "7 5 X2) (7) 

As long as the splitting is large enough, it is possible to 
have the decay X 2 —> Z'x l- For example, if the scalar 
giving rise to the Majorana mass is also the scalar re¬ 
sponsible for U{1)' breaking, M m can easily be of order 
M Z i. Here we have assumed a charge conjugation sym¬ 
metry, such that there is only one Majorana mass; if there 
are different Majorana masses for left- and right-handed 
components, diagonal couplings are also present. 

As in the previous model, we allow the Z' couplings to 
quarks and leptons to be set by additional free parame¬ 
ters. Our assumption is that the \2 has 100% branching 
to X\Z' , an( I that the Z' has 100% branching to qq , giving 
the final state signature of a dijet resonance plus miss¬ 
ing transverse momentum. For the dilepton plus missing 
transverse momentum signature, we allow for a signifi¬ 
cant branching fraction of the Z' to muons. 

To avoid scanning over too many parameters, we con¬ 
sider two sets of benchmarks for Mi. 2 . Since the cross 
section increases with lower %i mass, we include one op¬ 
timistic case with very light xi : 

Mi = 5 GeV, M 2 = Mi + M z > + A; A = 25 GeV (8) 

This case is somewhat tuned for large Z' mass, since 
it requires a cancellation between Dirac and Majorana 
masses. 

We also include a case where the fermion masses scale 
with M z >: 


Mi = M z ,/ 2, M 2 = 2 M z , (9) 

With Mi < M Z i , the interactions above are not sufficient 
for xi to obtain the correct thermal abundance in the 
standard cosmology. Since this is model-dependent, we 
leave this an open question and instead focus here on 
lighter dark sector masses, where the LHC sensitivity is 
better. 


Light Z' with Inelastic EFT coupling 

The models thus far rely on the Z' coupling to quarks 
in order to be produced at the LHC. Rather than produc¬ 
ing dark sector states through the new Z', we consider 
the possibility that it is produced through a new contact 
interaction: 

^2 ?7^ (x'27 / "7 5 Xi + Xi7 m 7 5 X2) • (10) 

Similar to the model just discussed, we have assumed 
two dark sector states Xi ,2 with an off-diagonal coupling 
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darkHiggs darkHiggs 




Light Vector 



Inelastic Vector 



Light Vector 



FIG. 3: Distribution of reconstructed (left) and rrijj (right) in the jj +J?t final state for each of the three models considered. 
We show a subset of our Z' mass points and consider the two cases for the masses of the other states, as discussed in the text. 
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to the Z'. The Z' + process is analogous to that 
of the previous section; however, we have effectively re¬ 
placed the intermediate s-channel Z' with a heavy Z' Hl 
where the Z' H has been integrated out to give the oper¬ 
ator above. For our benchmarks, the mass spectrum of 
the states is taken be the same as in Eqs. [8] and [9] 

The Z' produced in the decay can then be very weakly 
coupled to SM fermions, evading many direct search con¬ 
straints. For example, this small coupling could be gen¬ 
erated by kinetic mixing of the Z' with hypercharge and 
kinetic mixing parameter e 1. The only requirement 
is that the Z' decays to the visible fermions on collider 
time scales, which is easily satisfied for e > 10 -5 . For 
each search channel we show results assuming either a 
100 % branching fraction to jj or pp in order to match 
our signal regions; however the results can easily be 
scaled for the case of kinetic mixing where, for example, 
Br (pp) ss 0.12 for large Mz'- 

Similar ideas have been considered in hidden valley 
models Esna, which can give lepton jet signals from 
multiple light Z's [S3J [B4]. The main difference here is a 
looser signal requirement of a single Z' in the final state, 
and a wider range of Z' masses considered, which allow 
reconstruction of the dijet or dilepton resonance. 


LHC SENSITIVITY 

In the following sections, we consider the Z l —> jj and 
Z' —r It decay modes, propose an event selection and 
describe the expected sensitivity of the LHC dataset to 
Z' + for each of the models above. 


Dijet Mode 

Decays of a Z' to a pair of quarks results in two high- 
Px jets. In the following, the basic preselection requires 
at least two jets, each with px > 20 GeV and |? 7 | < 2.5. 
Events with a reconstructed electron or muon with px > 
10 and |p| < 2.5 are vetoed. 

The candidate Z' is reconstructed from the leading two 
Pt jets. To suppress the non-peaking backgrounds, a 
mass window rrijj £ [0.8 x mz' , mz' + 30 GeV] is applied. 
Distributions of rrijj and ^?x for the signal are shown in 
Fig-i For further details on how these distributions vary 
among the models, see the Discussion section. 

The primary background processes are Z —► vv in asso¬ 
ciation with two initial-state jets, or W —> tv in associa¬ 
tion with two initial-state jets and where the charged lep¬ 
ton is not identified. Events are simulated at parton level 
with MADGRAPH5 [55], with PYTHIA [BB] for showering 
and hadronization and DELPHES |S7j with the ATLAS- 
style configuration for detector simulation. Backgrounds 
are normalized to leading-order cross sections; the un¬ 
certainty is calculated by varying the factorization and 



FIG. 4: Distribution of reconstructed rrijj in the jj + 
final state, for the expected SM background as well as several 
examples of the signal in the dark Higgs (DH) model. Events 
are required to satisfy the preselection as well as have > 
300 GeV and leading jet px > 250 GeV, but no rrijj selection 
is applied. 


renormalization scales by factors of two. We validate 
our background model by comparing to the ATLAS re¬ 
sults [5] with rrijj £ [50,120] GeV and $x > 350,500 
GeV; the comparison is not precise due to the differences 
in the jet algorithm and radius parameters, but the esti¬ 
mates are roughly consistent. In Fig. [4] distributions of 
rrijj are shown for the expected backgrounds. 

To suppress the large dijet background, large is re¬ 
quired. The value of the threshold in |?x is determined 
by optimizing with respect to the expected upper lim¬ 
its on the cross section. In the case of the dark Higgs 
and light vector models, which have similar ^?x distri¬ 
butions, the threshold is ^?x > 200(300) GeV for values 
of mz' < 100 (> 100) GeV. In the case of the inelas¬ 
tic EFT model, which has larger the threshold is 
> 300(400) GeV for values of m z • < 100 (> 100) 
GeV. In addition, we require the px of the leading jet to 
be at least (^ x thresh — 50) GeV, which helps in suppress¬ 
ing the U+jets background. The efficiency of the final 
selection is shown in Fig [5] and detailed in Table [I] for 
various Z' and dark matter masses. 

Upper limits are calculated in counting experiments, 
using a profile likelihood ratio JBS] with the CLs tech¬ 
nique mm- Limits on the production cross section of 
cr(pp -> Z'xx -t jj'XX) are shown in Fig. [bj 


Dilepton mode 

Leptonic decays of a Z' may result in two high-px elec¬ 
trons or muons. In the following, the basic preselection 
requires at least two opposite-sign electrons or muons, 
each with px > 30 GeV and \rj\ < 2.5 as well as |?x >100 
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FIG. 5: {Top) Efficiency of the jj + selection described 
in the text, for two choices of mass spectra in each of the 
three models considered. Note that the minimum required 
increases above m z > > 100 GeV. ( Bottom ,) 95% CL upper 
limits on the production of Z' —¥ jj + as a function of the 
Z' mass. 

GeV and pt{££) > 80 GeV. Events with a third charged 
lepton or at least one jet with pt > 20 GeV and |?y| < 2.5 
are vetoed. 

Due to the tight constraints on Z' coupling to elec¬ 
trons discussed above, we will focus on the muonic chan¬ 
nel here. To a good approximation, the backgrounds 
would be larger by a factor of 2 if both lepton final states 
were included, and for models where the Z' decays to 
both charged lepton flavors, the resulting limits would 
be stronger by up to a factor of a /2 if systematic uncer¬ 
tainties are not dominant. 

The candidate Z' is reconstructed from the two lep¬ 
tons. To suppress backgrounds which do not peak at the 
Z' mass, a requirement that mu £ [0.9 x mz',mz' + 
25 GeV] is applied. Distributions of mu and are 
shown in Fig. [6j the dependence of these on different 
models and mass parameter choices is examined further 
in the Discussion section. 

The primary background processes are diboson pro¬ 
duction, such as ZZ —f ££vv, WZ —> £vt£, WW —> £v£v 


TABLE I: Signal efficiency and expected background yields 
for several Z' masses in the jj + final state. Only the 
heavy mass spectrum choice is listed. The background un¬ 
certainty is 27% obtained by varying the renormalization and 
factorization scales by factors of two. 



m z > [GeV] 

50 200 400 

[GeV] 

> 200 > 300 > 300 

Dark Higgs 

Light Vector 

Signal Efficiencies 
0.01 0.02 0.10 

0.002 0.03 0.20 

Z -A vv + jj 

W —>■ £u + jj 

Total Background 

Background Estimates 
3000 2,200 2,000 

350 300 330 

3,350 2,500 2,300 

[GeV] 

> 300 > 400 > 400 

Inelastic EFT 

Signal Efficiencies 
0.007 0.07 0.16 

Z -A vv + jj 

W —> tv + jj 

Total Background 

Background Estimates 
60 360 470 

10 50 65 

70 410 535 


or Z 7 — > t£vv. Top pair backgrounds are effectively sup¬ 
pressed via the jet veto. Events are simulated at par- 
ton level with MADGRAPH5 [63], with PYTHIA f66] for 
showering and hadronization and DELPHES [67] for detec¬ 
tor simulation. Backgrounds are normalized to leading- 
order cross sections; the uncertainty is calculated by 
varying the factorization and renormalization scales by 
factors of two. A minimum 15% uncertainty is applied 
to cover uncertainty due to the high-px region. We vali¬ 
date our background model by comparing to the ATLAS 
results [13[| with mu £ [76,106] GeV and 1/Ht > 150, 250, 
350, 450 GeV; our estimates agree within uncertainties. 
In Fig. [TJ distributions of mu are shown with the ex¬ 
pected background. 

As in the dijet case, large missing transverse momen¬ 
tum is required to suppress the large ££ backgrounds; the 
requirement ^?x > 100 is found to give the strongest ex¬ 
pected limits across all models and masses. The efficiency 
of the selection is shown in Fig [8] and detailed in Table |TT| 
for various Z' and dark matter masses. 

Upper limits are calculated in counting experiments, 
using a profile likelihood ratio [63] with the CLs tech¬ 
nique uni ED]- Limits on the production cross section of 
a{pp — > Z'xx P + XX) are shown in Fig. [8j 

DISCUSSION 

The kinematic distributions in |?x and invariant 
masses of the different models are shown in Figs. [3] and [6] 
In both the dark Higgs and light vector models, the inter- 

















Fraction of Accepted Events Fraction of Accepted Events ^ Fraction of Accepted Events 


darkHiggs 


darkHiggs 


Z' Mass (GeV) 



Light Vector 




Light Vector 



Inelastic Vector 
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FIG. 6: Distribution of reconstructed (left) and mu (right) in the ££+$t final state for each of the three models considered. 
We show a subset of our Z' mass points and consider the two cases for the masses of the other states, as discussed in the text. 
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FIG. 7: Distribution of reconstructed met in the p + p~ +$t 
final state, for the expected background as well as a signal 
example. The IFT label refers to the inelastic EFT model. 
Events are required to satisfy the preselection as well as have 
> 100 GeV and pee > 80 GeV, but no met selection is 
applied. 


TABLE II: Signal efficiency and expected background yields 
for several Z' masses in the p + pT +|?t final state with > 
100 GeV. In each model, the masses are are chosen to be that 
of the heavy spectrum case. 



m z i [GeV] 

50 200 400 

Model 

Signal Efficiencies 

Dark Higgs 

0.06 

0.13 

0.17 

Light Vector 

0.01 

0.14 

0.18 

Inelastic EFT 

0.09 

0.16 

0.18 

Process 

Background Estimates 
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mediate s—channel Z' is off-shell, and so the spectra 
are typically softer than in the inelastic EFT model and 
primarily determined by total mass in the final state. As 
a result, the high |?x tail can look similar for different Z' 
masses, if the other masses are correspondingly adjusted. 
Note that for the dark Higgs model, the spectra de¬ 
pends on the mass of the dark Higgs and Z ', and not 
directly on the dark matter mass, while for the light vec¬ 
tor model the spectra depend on the total mass in the 
X1X2 final state as well as on their splitting. 

In the inelastic EFT model, production goes through 
a higher dimension operator, leading to harder |?x spec¬ 
tra and less sensitivity to the masses in the final state. 
Note that the high |?x tail in the Mz> = 50 GeV case 




FIG. 8: ( Top) Efficiency of the p + p~ +$t selection described 
in the text as a function of the Z' mass, for two choices of 
mass spectra in each of the three models considered. ( Bottom ) 
95% CL upper limits on the production of ( Z' —> p + p~) +$t 
as a function of the Z' mass. 


has an additional suppression, however, since such highly 
boosted low-mass Z' are unlikely to be resolved as two 
individual jets. Another effect that becomes important 
is the size of the splitting m X2 — m Xl compared to Mz>: 
when the xi is very light and the splitting is very to close 
to Mz' , the px of \2 is transferred nearly entirely to the 
Z' and consequently the ^?x spectrum is harder. This 
corresponds to the case in Eq. [8j Conversely, less px is 
transfered to the Z' as the splitting is increased and as xi 
becomes heavier, as in Eq. [9] The effect competes against 
the increase in missing transverse momentum with larger 
m X2 , m Xl . For the cases shown here, as we increase the 
dark matter masses, we also scale the splitting up accord¬ 
ingly. As a result, for a given Z ', the ^x distribution does 
not change much for the two different mass spectra we 
consider. 
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FIG. 9: Expected upper limits at 95% CL on the product of couplings g q gh D as a function of M z / for the dark Higgs model, for 
8 TeV pp collisions in two different mass benchmarks. Left, the sensitivity of the jj + $t channel is compared to the constraint 
on the hadronic Z width (labelled Rz), shown in black for gh D = 1 (solid) and gh D = 2 (dashed), as well as direct dijet 
resonance searches m for a new Z'. Right, the sensitivity of the pp + t channel is compared to various dimuon resonance 
searches at CDF 33 and ATLAS [ 50 ], all shown for gh D = 1. The low-mass dimuon limits are interpreted from the results of 
Ref. |54 |: both 7 TeV recast limits (dotted) and 8 TeV sensitivity projections (dashed) are shown. We do not consider masses 
in the grey shaded region due to the extremely large Drell-Yan background near the Z mass. 



50 


100 


800 


M z , [GeV] 



M z , [GeV] 


FIG. 10: Expected upper limits at 95% CL on the product of couplings g q gh D as a function of M Z ' in the Light Vector model, 
for 8 TeV pp collisions in two different mass benchmarks. The dijet and dilepton resonance limits are the same as those in 
Fig. § with g x — 1 for all dilepton resonance limts. 


Model Constraints 

We evaluate the sensitivity of the first LHC run to 
each of the models presented in this paper. The results 
are shown in Figs. m considering both dijet and dilep¬ 
ton resonances in the mass range Mz’ = 50 — 800 GeV. 
For each final state, we show results assuming a 100% 
branching ratio of the Z' to dijets or to dimuons accord¬ 
ing to our signal regions. For the dimuon final states, 
we do not consider the mass range Mz' £ (65,100) GeV 
since there is a significant Drell-Yan background at these 


invariant masses, as shown in Fig. [7} 

As discussed in the Constraints section, there are 
strong constraints on electron couplings to the Z ', which 
severely limits the Z' branching ratio to electrons in the 
dark Higgs and light vector models. For uniformity in 
our presentation of results we have therefore considered 
only the muonic final state. The combined dimuon and 
di-electron result would be somewhat stronger in the case 
that the Z' decays to both flavors equally, as in the in¬ 
elastic EFT model. 

Constraints for the dark Higgs model are shown in 
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FIG. 11: Expected lower bound at 95% CL on A from in the Inelastic EFT model, for 8 TeV pp collisions in two different 
mass benchmarks. The branching ratio of the Z' to jets and muons is taken to be 100% in each case. 


Fig- i for each of the two choices of dark Higgs mass 
given in Eqs. |4][5j The predicted cross section for the 
Z' +|?t signal is proportional to g\ D g\- For M z > < 200 
GeV, the constraints for the two Mh D cases are similar. 
The lighter Mh D case has a larger cross section, but at 
the cost of a softer spectrum and hence reduced se¬ 
lection efficiency, as shown in Fig. [5] and Fig. [8] Above 
M Z ' = 200 GeV, the limits on the heavy Mh D scenario 
become significantly weaker due to the rapidly decreasing 
production cross sections. 

The missing transverse momentum searches are com¬ 
pared in each case with the corresponding direct dijet or 
dilepton resonance searches from various hadron collid¬ 
ers. Since the Z' + limits depend on an additional 
model parameter gh D , we show the resonance search lim¬ 
its for a reference value of gh D = 1; if this coupling were 
stronger, these limits would be relatively weaker. As can 
be seen, for M Z ' > 150(100) GeV in the dijet (dilep¬ 
ton) case, the direct resonance searches give stronger con¬ 
straints on the model. 

At low M Z f , constraints from the experimental col¬ 
laborations are not available. However, we compare 
the dimuon results with the low mass dimuon resonance 
study in Ref. t 54] , finding that their recast limits of 7 TeV 
data would still be stronger than that from Z' +$t- Al¬ 
though a search helps reduce backgrounds, the statis¬ 
tics for the signal are also lower: in this model the mono- 
Z' signal requires an off-shell intermediate Z' and the 
production of an additional particle (the dark Higgs) in 
association with the Z ', thus leading to a suppression of 
~ 10 3 in the rate even for the “light” M hD case. 

We find the most relevance for this signal model in 
the context of leptophobic Z's with mass below ~ 150 
GeV, where there is a gap in existing dijet resonance 
studies. As discussed in the Constraints section, there is 
an indirect constraint since a light Z' would modify the 


hadronic Z width, which we show in Fig. [9] for gh D = 1 
and gh D = 2. An LHC associated Z' search ""[28] offers the 
best prospects for robust constraints competitive with 
the Z' + results in this mass range. 

The limits in the light vector model are shown in 
Fig. 10 and the behavior is qualitatively similar. In 
addition, we make the analogous assumptions as in the 
dark Higgs results described above, with the resonance 
search results shown for g x = 1. We find that the dijet 
resonance plus performs more favorably here, hav¬ 
ing the best sensitivity to the light m Xl scenario below 
M z > ss 200 — 300 GeV. However, the dimuon plus miss¬ 
ing transverse momentum search would again be weaker 
than a direct dimuon search in the entire mass range. 

Finally, the inelastic EFT model limits are shown in 
Fig- EH We constrain A, the scale of the operator leading 
to dark matter production, for each of the two channels. 
Since the Z' can be very weakly coupled in this model, 
the dijet and dimuon resonance limits above do not apply 
and by construction, the Z' +$t search provides the best 
constraint. This model is especially interesting for the 
dimuon mode, where limits on A reach roughly 5 TeV, or 
around 3 TeV if rescaled to Br(/r + /i _ ) = 0.12. 

We also compare the results of our Z' plus missing 
transverse momentum search to constraints derived us¬ 
ing existing ^x-based searches for new physics beyond 
the standard model. For the dijet resonance plus missing 
transverse momentum case, the monojet search region 
would be sensitive to our models since up to two jets are 
allowed. However, by focusing on specific m 7J - windows, 
our analysis has far lower backgrounds. We also com¬ 
pare with the multijet plus missing transverse momen¬ 
tum SUSY search jJT]: we find the SUSY study is less 
sensitive to our models, since it requires a much larger 
amount of visible and missing transverse momentum in 
order to optimize for a signal from new heavy colored par- 
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tides. In the dilepton resonance case, we compare with 
the cliargino search m- Here we find a fair amount of 
overlap in the signal regions, leading to comparable sen¬ 
sitivity to our models; for a more detailed discussion of 
the bounds obtained from applying the chargino search, 
see [361 |. 


CONCLUSIONS 

We have presented a new collider signal for dark mat¬ 
ter: missing transverse momentum and a dijet or dilep¬ 
ton resonance. This work adds to the existing mono-X 
and simplified models of missing transverse momentum 
signals, expanding the coverage of LHC searches to new 
dark sector physics that may be difficult to observe in 
other channels. In this paper, we introduce several sim¬ 
plified models for a Z' produced in association with the 
dark matter, determine the sensitivity of the current LHC 
dataset to these models, and compare with other collider 
searches for Z' s. 

When the Z' plus dark matter production relies on the 
Z' couplings to quarks, we find that a mono -Z' channel is 
more sensitive than dijet resonance searches only below 
Mz' of a few hundred GeV. In this mass range, there are 
currently no published results searching for a resonance 
from a hadronically decaying Z ', and the requirement of 
can significantly reduce the QCD dijet background. 
On the other hand, in these same models, when the Z' 
can also decay to leptons then a direct dilepton resonance 
search is expected to be a more powerful constraint on 
the model in the entire Z' mass range. 

The Z' can also be produced in the decay of dark sector 
states, which are coupled to quarks through an effective 
contact interaction. Then the Z' may be weakly coupled 
to SM states, easily satisfying other direct collider con¬ 
straints. Such a model would be challenging to observe 
in other missing transverse momentum searches, but give 
rise to a mono -Z' signal. As the first run of the LHC has 
shown, there is need for a broad range of dark matter sig¬ 
nals to explore the many possibilities for the dark sector 
and to take full advantage of the data. 
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